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We find a new Tully-Fisher-like relation for spiral galaxies holding at dif- 
ferent galactocentric radii. This Radial Tully-Fisher (RTF) relation allows us 



. to investigate the distribution of matter in the optical regions of spiral galax- 

ies. This relation, applied to three different samples of rotation curves of spiral 

galaxies, directly proves that: 1) the rotation velocity of spirals is a good mea- 
*sO '. 

, sure of their gravitational potential and both the RC's amplitudes and profiles 

. are well predicted by galaxy luminosity 2) the existence of a dark component, 

less concentrated than the luminous one, and 3) a scaling law, according to 

o ■ 

?H ' which, inside the disk optical size: Md ar k/Mi um = 0.5(Lb/10 11 Lb ( t v ,; ' 

C3 

•Si 1 INTRODUCTION 

In 1977 Tully and Fisher discovered that the maximal rotational velocity V max of a spiral 
galaxy, measured by the FWHM of the neutral hydrogen 21-cm line profile, correlates with 
the galaxy luminosity by means of a power law of exponent a ~ 4. This equivalently reads 
as: 

M = ahgV max + b (1) 

with M the absolute magnitude in some specified band and b a constant. It was immediately 
realized that this relation, hereafter TF, could serve as a powerful tool to determine the dis- 
tances of galaxies (Pierce & Tully 1988) and to study their dynamics (Persic & Salucci 1988). 
The rotational velocity reflects the equilibrium configuration of the underlying galaxy grav- 
itational potential, especially when V max is directly derived from extended rotation curves. 
Before proceeding further, let us point out that spiral galaxies have a characteristic size 
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scale, R op t, that sets also a characteristic reference velocity V(R opt ). Rd, the exponential 
thin disk length scale, is a natural choice for such reference radius; in this paper, however, 
we adopt for the latter a minimal variant, i.e. a multiple of this quantity: R opt = 3.2Rd 
(see Persic, Salucci & Stel 1996, hereafter PSS). (No result here depends on the value of the 
multiplicity constant). This choice is motivated by the fact that 3.2R D , by enclosing 83% 
of the total light, is a good measure of the "physical size" of the stellar disk, and that, for 
many purposes, V opt = V(R opt ). * 

Let us stress that some known kinematical quantities are not suitable reference veloci- 
ties. For example, the value of V max for a spiral depends on the extension and on the spatial 
resolution of the available RC and, in addition, it does not have a clear physical interpre- 
tation, sometimes coinciding with the outermost available velocity measure, in other cases 
with the innermost one. Also Vi ast , the velocity at the outermost measured point obviously 
does not have a proper physical meaning, in addition some spirals never reach the, so called, 
asymptotic flat regime (PSS and Salucci & Gentile 2006). 

Coming back to the TF relation, its physical explanation, still not fully understood, 
very likely involves the argument that in self-gravitating rotating disks both the rotation 
velocity and the total luminosity are a measure of the same gravitational mass (e.g. Strauss & 
Willick 1995). Notice that, if this argument is correct, both V max and Vi ast are just empirical 
quantities of different and not immediate physical meaning. 

The stars in spiral galaxies are settled thin disks with an exponential surface mass dis- 
tribution 

E(R) = Z e- R/RD , S = hL s , (2) 

where S = (M d /L)I is the central surface mass density, with I the central surface bright- 
ness, in the first approximation, constant in spirals. L is the total luminosity in a specific 
band, k\ and s are constants. Observationally, as a good approximation, we find that (e.g. 
PSS): 

Rd = k 3 L q , (3) 

with q = 1/2, and k^ = const. It is illustrative to set, for the time being, s = 0, i.e. to assume 
that the disk mass-to-light ratio has an universal value. Let us consider the condition of self- 
gravity equilibrium for the stellar disk, i.e. the ratio E = GM D /{V^ pt R opt ). By writing: 

* For the PS95 sample: log V opt = (0.08 ± 0.01) + (0.97 ± 0.006) \ogV max - 
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E = k 2 l} , 



(4) 



where k 2 and t are constants, we have that Freeman disks are completely self-gravitating 
for k 2 — 1-1, t — 0. Combination of the previous equations and the above assumptions leads 
to the well known relation: luminosity oc (velocity) 4 . Random departures from the above 
conditions will produce a larger scatter in the TF relation, while systematic departures from 
the assumption made above e.g. variations of the stellar population with galaxy luminosity 
or violation of the condition of self-gravity, will modify the slope, zero-point and scatter, 
possibly in a band-dependent way. In fact, by relaxing some of the assumptions made above, 
we have the more general relationship: L oc V^} +s+q ^ (here s, q and t can be band- 
dependent) that can be even more complex and non-linear when the scaling laws (2), (3), 
(4) are not just power laws. As a matter of fact, in several different large samples of galaxies 
it has been found that the TF has different slope and scatter in different bands: ai ~ 10, 
si ~ 0.4mag, while a B ~ 7.7, s B ~ 0.5mag (Pierce & Tully 1992), (Salucci et al. 1993). 
Moreover, a non linearity in the TF is often found at low rotation velocities (Aaronson et 



We know that spiral galaxies are disks of stars embedded in (almost) spherical halos 
of dark matter and this is crucial for understanding the physical origin of the TF relation 
(Persic & Salucci 1988, Strauss & Willick 1995, Rhee 1996). It is well known that the dark 
halos paradigm is supported by the (complex) mass modelling of galactic rotation curves 
(e.g. PSS and references therein) and it implies that disks are not fully self-gravitating. At 
any radius, both the dark and luminous components contribute to the (observed) rotational 
velocity V(R), with a relative weight that varies both radially and from galaxy to galaxy. The 
resulting model circular velocity can be written as a function of the useful radial coordinate 
x = Rj Ropt as t : 



al. 1982). 



V model (x) = (GM D /R D y/ 2 [f d {x)/f d {l) + Tf h {x,a)Yl\ 



(5) 



where fd(x) is the Freeman velocity disk profile 




(6) 



and 



f h (x,a) = ( 



x 



2 + a 



^)(l+« 2 ) ? 



(7) 



t For simplicity here we neglect the bulge, we will consider it in section 4. 
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Figure 1. The TF relation at R op t (triangles); the predicted relation at « (l/3)R op t in the case of: no-DM (short dashed 
line), T = 1 NFW halo (dotted line) and a = 2 and T = O.5(L S /1O 11 L S0 ) (solid line). 

M d is the disk mass, T is the dark/visible matter velocity ratio at R opt and a is the halo 
velocity core radius in units of R opt . The adopted halo function f h (x,a) (see PSS) is the 
simplest way to describe the contribution of dark matter halo; in fact, for an appropriate 
value of a the DM term in V moo iei(x) describes both the "empirical" universal rotation curve 
halo velocity profile in PSS and, in the regions under study (by setting a = 1/3) the 
"theoretical" AC DM Vnfw(x) Navarro- Frenk- White halo velocity profile. 
Once a spiral mass model is assumed, the well know TF relation at R op t- 



(e.g. PSS, Persic & Salucci 1995 hereafter PS95), may lead through eq.(5)-(7), to the ex- 
istence of a family of similar relationships at radii jR op t {j — 0.2 — 1.4), whose values for 
the slopes, scatters and zero points are all different and all depending on the characteristics 
of the mass distribution. In Fig.l we show this for the reasonable values of T and a at. Of 
course if T and a are random functions then no relation will emerge. 

The aim of this paper is to investigate the rotation curves (RCs elsewhere) of three 
samples of spiral galaxies and to extract from actual data the Radial Tully-Fisher relation 
(RTF), i.e. a family of TF-like relations holding at any properly chosen radial distance. 




(8) 
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Then, we will use the properties of such relationships to investigate the mass distribution in 
these objects. 

Let us notice that the strict Universal Rotation Curve paradigm, put forward in PSS, 
implies, in addition to an Universal velocity profile at any chosen luminosity, also a "TF 
like" relationship at any chosen radius x. While in PSS the aim was mainly to derive the 
above universal profile, in the present paper, meant to be complementary to PSS, we will 
head for the latter set of relationships. 

In order to obtain statistically relevant results as free as possible from biases, observa- 
tional errors and cosmic variance effects, we use three different samples of spiral galaxies with 
available rotation curves. These samples mostly contain Sb-Sc spiral galaxies. The number 
of early type spirals and dwarfs is very small. Moreover, in general, the bulge affects only the 
first reference radius. Different investigation will be necessary to assess the present results 
in bulge dominated spirals and HI dominated dwarfs. 

The plan of this work is the following. In section two we describe our data samples 
and the main steps of the analysis. In section three we show that TF-like relations hold 
at specific radii, and we derive the basic parameters of these relations for our samples. In 
the next section we discuss the implications of the existence of Radial TF relation and we 
propose a simple mass model that fits the data. The conclusions are given in section 5. 

2 DATA AND ANALYSIS 

Sample 1 consists of 794 original RCs of PS95 (notice that for most of them the limited 
number of independent measurements and sometimes some circular motion make difficult 
to derive a proper mass model, but instead will be possible with the method in section 2,3). 

In each RC the data are radially binned on a 0.2R opt scale so that we have 4-7 independent 
and reliable estimates of the circular velocity, according to its extension. 

Sample 2 from Courteau (1997) consists of 86 RCs (selected from 304 galaxies) and 
Sample 3 of Vogt (2004) 81 RCs (selected from 329 galaxies). These samples have been built 
by selecting from the original samples only objects with high quality and high resolution 
kinematics yielding reliable determinations of both amplitudes and profiles of the RCs. To 
ensure this, we have set the following selection criteria. The RCs must: (a) extend out to 
~ R op t', (b) have at least 30 velocity measurements distributed homogeneously with radius 
and between the two arms; and (c) show no global asymmetries or significant non-circular 
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motions: the profiles of the approaching and receding arms must not disagree systematically 
more than 15% over lR d length-scale. The velocity errors are between 1% — 3%. 

In each galaxy we measure the distance from its center R in units of R opt (R op t = S.2Rd) 
and we consider a number of radial bins centered at R n = (n/5)R opt for the PS95 sample 
and at R n = (n/20)R opt for the other two samples; we take the bin size 5 = 0.2R opt for the 
PS95 sample and 5 = 0.0QR opt for the other two samples. Then we co-add and average the 
velocity values that fall in the bins, i.e. in the radial ranges R n — 5/2 < R n < R n + 5/2 
and we get the average circular velocity V n at the chosen reference radii R n . (For the PS95 
sample this was made in the original paper) 

A "large" radial bin size has been chosen for the (much bigger) PS95 sample because, 
by selection, most of its RCs have a relatively smaller number of measurements. For the 
other two samples, that, exclusively include extended high quality RC and large number of 
measurements we decrease the bin size by a factor 3.2. 

In short, we will use two different kind of samples: sample 1 includes 794 Sb-Sd galaxies 
with I magnitudes whose RCs are estimated inside large radial bins that smooth out non 
circular motions and observational errors; samples 2 and 3 include 167 galaxies with R 
magnitudes, whose RCs of higher quality are estimated inside smaller radial bins providing 
so a larger number of independent data per object. 

We look for a series of correlations, at the radii i?„ between the absolute magnitude M 
(in bands indicated below) and log V n = log V(R n ). Data in the / (Mathewson et al. 1992) 
and r (Courteau 1996), (Vogt et al. 2004) bands will allows us to check the dependence of 
our results on the type of stellar populations in spiral galaxies. Finally let us stress that the 
uncertainties of photometry are about 10% and therefore negligible. 



3 THE RADIAL TF RELATIONSHIP 

Given a sample of galaxies of magnitude M band and reliable rotational curves, the Radial 
Tully-Fisher (RTF further on) relation is defined, as the ensemble of the fitting relationships: 

M hand = a n log V n + b n , (9) 

with a n , b n , the parameters of the fits and R n the radial coordinates at which the relationship 
is searched. The latter is defined for all object as a fixed multiple of the disk length-scale (or 
equivalently of R op t)- Parameters a n , b n are estimated by the least squares method (without 
considering the velocity /magnitude uncertainties). 
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Figure 2. The Radial Tully-Fisher relations for the PS95 sample. Each one of the 6 relations is indicated with different color. 

The existence of the Radial TF relation is clearly seen in Fig. 2 and Fig. 3, where all the 
TF-like relations for the PS95 sample are plotted together and identified with a different 
color. It is immediate to realize that they mark an ensemble of linear relations whose slopes 
and zero-points vary continuously with the reference radius R n . 

Independent Tully-Fisher like relationships exist in spirals at any "normalized" radius 
R n . We confirm this in a very detailed and quantitative way in the Figures 11, 12 and in the 
Tables 1, 2, 3, where very similar results are found for the other two samples. It is noticeable 
that the various investigations lead to the same consistent picture. 

The slope a n increases monotonically with R n ; the scatter s n has a minimum at about 
two disk length-scales, 0.6R opt . In the J-band the values of the slopes are about 15% larger 
than those in the r-band. This difference, well known also for the standard TF, can be 
interpreted as due to the decrease, from the r to the I band, of the parameter s (see Eq. 
2), as effect of a different importance in the luminosity of the population of recently formed 
stars (Strauss & Willick 1995). 

It is possible to compare the RTFs in different bands; in the case of absence of DM, true 
in the inner regions of spirals (except in the very luminous galaxies and LSB almost absent in 
our sample), for the reasonable values s r = 0.1 and Sj = the power law coefficient of the Lj 
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Figure 3. The Radial Tully-Fishcr relation for the PS95 sample. 

vs velocity relationship is larger by a factor (0.5 + s r ) / (0.5 + s/) than that of the L r vs velocity 
relationship, in details by a factor 1.2. This correction allows us to compare the a n slopes as 
a function of R n for our samples (see Figure 4). Remarkably we find that the values of the 
slopes vary with R n according to a specific pattern: b n = —2.3—9.9(R n /R op t)+3.9(R n /R < y P t) 2 . 

It is worth to look at the scatter of the Radial TF relation (see Fig. 5). We find that, 
near the galactic center the scatter is large 0.3 — 0.4 dex, possibly due to a "random" 
bulge component governing the local kinematics in this region being almost independent of 
the total galaxy magnitude. The scatter starts to decrease with i? n , to reach a minimum of 
0.15—0.3 dex at R n corresponding to two disk length scales, the radius where the contribution 
of the disk to the circular velocity V(R) reaches the maximum. From 2R D onward, the scatter 
increases outward reaching 0.3 dex, at the farthest distances with available data, i.e. at 3-4 
disk length-scales. 

Let us notice that these scatters are remarkably small. Most of the relations in the RTF 
family are statistically at least as significant as the standard TF relation, while the most 
correlated relationship, i.e. that at IRd shows a r.m.s of only 0.2 - 0.3 magnitudes (besides, 
significantly smaller than that of the standard TF (see Figure 6)). Since the r.m.s include 
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Figure 4. The slope of the Radial Tully-Fisher relation at different radii for the 3 samples. The slope for the standard TF is 
about to -7.5. 
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Figure 5. The scatter of the Radial Tully-Fisher relation at different radii for the 3 samples. 

also the effects of various observational errors, such small values could even indicate that 
the intrinsic scatter of the RTF relation at ~ 2Ru is almost zero. 

An important consequence of the smallness of the r.m.s of TF-like relationships is that 
we can claim that at any radius, the luminosity predicts the rotation velocity within the 
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This error (that includes also distance, inclination, and non circular motions errors) is re- 
markably small, mostly lying between 5%-10%, and not exceeding 20% even in the very 
inner bulge dominated regions (see figure for the PS95 sample). This result is an additional 
proof for the URC paradigm, even more impressive than the set of synthetic RC's in PSS. 

Incidentally, the existence of a radius (R = O.QR opt ) at which the TF-like relations show 
a minimum in the internal scatter is not related to the overall capability of the luminosity in 
predicting the rotation velocities. At very small R ~ (0.1 — 0.2)R opt the (random) presence 
of a bulge increases the scatter, in that, (as in ellipticals) the actual kinematical-photometric 
fundamental plane likely includes a third observational quantity (maybe the central surface 
density). At large R, the observed small increase of the scatter can be avoided by introducing 
small quadratic term in the eq (9). 

The scatter of the RTF in the R band for the corresponding samples is somewhat larger 
than that in the I band for the PS95 sample. This can be easily explained by the following: 
i) the former samples include also a (small) fraction of Sa objects and their RCs are of 
higher spatial resolution (lower bin size) and therefore less efficient in smoothing the non- 
circular motion arisen by bars and spiral arms ii) the R band is more affected than the I 
Band by random recent episodes of star formation. A conservative estimate of these effects 
is cr obs > 0.2, thus the intrinsic scatter of the RTF in the R band (s 2 n — u^) - 5 results similar 
to that found in the I band. 

4 RADIAL TF: IMPLICATIONS 

The marked systematic increase of the slopes of the RTF relationship, as their reference 
radius R n increases from the galaxy center to the stellar disk edge, bears very important 
consequences. First, it excludes, as viable mass models, those in which: 

i) The gravitating mass follows the light distribution, due to a total absence of non 
baryonic dark matter or to DM being distributed similarly to the stellar matter. In both 
cases, in fact, we do not expect to find any variation of the slopes a n and very trivial 
variations of the zero points b n with the reference radius R n , contrary to the evidence in 
Tables 1, 2, 3 and Fig. 4. 

ii) The DM is present but with a luminosity independent fractional amount inside the 
optical radius. In this case, in fact, the value of the circular velocity at any reference radius 
R n will be a luminosity independent fraction of the value at any other reference radius 
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(i.e. logV(R n ) = k nm + logV(i? m ), k nm independent of luminosity). As a consequence the 
slopes a n in the RTF will be independent of R n , while the zero-points b n will change in a 
characteristic way. It is clear that the evidence in Table 1, 2, 3 and Figure 4 contradicts this 
possibility. 

The RTF contains crucial information on the mass distribution in spirals. In paper II we 
will fully recover and test it with theoretical scenarios. Here, instead, we will use a simple 
mass model (SMM), that includes a bulge, a disk and a halo mass component and it is tunable 
by means of 4 free parameters. By matching this model with the slopes of the RTF relation 
vs reference radii relationship (hereafter SRTF) we derive the gross features of the mass 
distribution in spirals. Let us point out, that this method has a clear advantage with respect 
to the mass modelling based on RCs. In this latter procedure, since the circular velocities have 
a quite limited variation with radius, physically different mass distributions may reproduce 
the observations equally well. Here instead, on one side, we will use an observational quantity, 
the slope of the RTF relation, that shows large variation with reference radius; on the other 
side, physically different mass distributions predict very different slope vs reference radius 
relationships. 

First, without any loss of generality affecting our results, we assume the well known 
relationships among the crucial structural properties of spirals: i) 

R D = R 1 l 05 , (Ha) 

see PSS, with I = io(Mi-m™*)/2-5 anc i Mf ax = -23.5 the maximum magnitude of our 
sample, and ii) 

M D = M 1 l 1 - 3 , (116) 

(e.g. Shankar et al. 2006 and references therein). Notice that the constants R\ and Mi will 
play no role in the following. We will best fit the a n data, i.e. the SRTF relation shown in 
Fig (4) with the model function <2smm(x) derived from the SMM. In this way we will fix the 
free structural mass parameters. 

We describe in detail the adopted SMM, the circular velocity is a sum of three contri- 
butions generated by the bulge component, taken as a point mass situated in the center, a 
Freeman disk, and a dark halo. 

yJ„„M = gm d /r dIW) + + ^AM], 

where M ha [ is the halo mass inside R opt - It is useful to measure Vg MM in units of GMi/Rx, 
and to set it to be equal 1. 
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The disk component from Eqs. (6) and (10), (11) takes the form: 



V d 2 (x,l) = P*f d (x), (12) 

with V^(l,l) = 0.347. 

We set Mf, the bulge mass to be equal to a fraction (c&/(3.2 • 0.347)/ ' 5 of the disk mass, 
with Cb a free parameter of the SMM and the exponent 0.5 is suggested by the bulge-to disk 
vs total luminosity relation found for spirals. Then, we get 

V b 2 (x) = c b V d 2 (l,l)l°- 5 x- 1 . (13) 

The halo velocity contribution follows the profile of Eq. 7; at R opt it is set to be equal to 
Ch/ (3.2-0.347)/( fc,1 ~ a5 ) times the disk contribution Cb/ch and k h , and a are the free parameters 
of the SMM. 

Then, we can write: 

V* MM (x, a, I) = (cbl^/x + l°- 8 f d (x) + c h l^-°^f h (x, a)) , (14) 

where fd is given by Eq. (6), fh by Eq. (7). Notice that the simple form of Vsmm allows us 
to get the predicted slope parameter clsmm(x)- 

The core radius a is a DM free parameter, however let us anticipate that, provided that 
this quantity lies between 0.5 and 2 (see Donato et al. 2004) it does not affect in a relevant 
way the SMM predictions. Instead the a S MM(%) relationship strongly depends on the values 



© 0000 RAS, MNRAS 000, 000-000 



14 Yegorova & Salucci 





1.1 






1.0 


^ 1=3 

^^^^ 1=0.03 




0.9 






8 




a 

+ 


0.7 


* 




0.6 






0.5 






0.4 






0.3 






0.2 





00 02 OA 06 08 " L0 * L2~ 



Figure 9. Baryonic mass fraction as a function of normalized radius for high and low luminosity objects. 
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Figure 10. Baryonic mass fraction at two different radii as a function of luminosity. 

of Ch, kh, Cb, and therefore they can be estimated with a good precision. We can reproduce the 
observational a n = a(R n ) relationship, by means of the SMM (see Fig. 8) with the following 
best fit parameters values: i) kh = 0.79 ± 0.04, that means that less luminous galaxies have 
larger fraction of dark matter, ii) q, = 0.13 ± 0.03, and iii) Ch = 0.13 ± 0.06 that indicates 
that at (l,x) = (1, 1) (i.e. inside R op t), 20% of the mass is in the bulge component, 20% in 
the halo, while 60% is in the stellar disk. The quoted uncertainties are the formal \ 2 fitting 
uncertainties. 
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5 DISCUSSION AND CONCLUSIONS 

In spirals, at different galactocentric distances measured in units of disk length-scales jRp 
(j = 0.2, ...,4), there exists a family of independent Tully-Fisher-like relationships, M band = 
bj + cij log V(Rj), that we call the Radial Tully-Fisher relation, that contains crucial infor- 
mation on the mass distribution in these objects. In fact: 

1) The RTF relationships show large systematic variations in their slopes aj (between 
—4 and —8) and a r.m.s. scatter generally smaller than that of the standard TF. This rules 
out the case in which the light follows the gravitating mass, and, in particular, all mass 
models that imply: a) an absence of dark matter, b) a single mass component c) the same 
dark-to-luminous-mass fraction within jRo in all galaxies. 

2) The slopes aj decrease monotonically with Rj, implies the presence of a non luminous 
mass component whose dynamical importance, with respect to the stellar disk, increases 
with radius. 

3) The existence of the RTF and the features of the slope vs Rj can be well reproduced 
by means of a three components mass model that includes: a dark halo, with a core of radius 
4Rd and mass Mhaio(Ropt) oc I - 7 , a central bulge with Mb u i ge oc / L8 , an exponential thin disk 
of mass Ma oc I 13 with, at 10 11 Lb,©, 60% of the mass inside R opt in stellar form. 

Let us also notice that we have produced a qualitatively new evidence for the presence 
of a luminosity dependent mass discrepancy in spirals, different from that obtained from 
the non Keplerian shapes of the RC's. While, the latter originates from a failure: we do not 
observe the Keplerian fall-off of the circular velocities at the disk edge, and therefore we 
must postulate a new component, here, we provide a positive evidence for the existence 
of such dark component: we detect radial change of the slope and the scatter of existing 
relations between observables that positively indicates the presence of a more diffuse dark 
component. 

The small scatter of the TF-like relationships in the inner regions of the spirals (x = 0.4 — 
0.6) has important implications for the DM halo core vs cusp controversy, more specifically 
on the claim that the observed rotation velocities V(R) do not coincide with the circular 
velocities, i.e. with the equilibrium velocities associated with the central galaxy gravitational 
potential $, V C i rc (R) = (—RdQ/dR) 1 / 2 ). It is claimed that, due non-circular/non-central 
motions, gas pressure gradients and other effects, V(R) < V C i rc by a substantial amount 
and models in which observed cored-solid body rotation curves V(R) underly actual NFW- 
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cusped circular velocities V C i rc (R) = Vnfw{R) 7^ V(R) have been built (Hayashi 2006). Let 
us notice that this effect must be very large in that it must explain the large discrepancy 
between the observed rotation curves (taken as bona-fide circular velocities) and their NFW 
best fit mass models: the typical difference is 

(V - V NFW ) = (20 - 60)A;m/s 

for galaxies with V opt = 100 — 150km/s (Gentile et al., Donato et al. and references therein). 
At R n = (0.4 — 0.6)R opt we have s n ~ 0.3, that certainly includes a component due to 
distance and inclinations errors. Even if we assume that its all due to the fact that V ^ V C i rc 
this would trigger a discrepancy of only ~ 5 — 15 km/s, far too small to account for the 
severe discrepancies of the NFW models. 

Finally, let us stress that any model of formation of spiral galaxies must be able to 
produce (e.g. in the I band) a Mj vs log V(2Rd) relationship with a slope of 7 ± 0.1 and an 
intrinsic scatter of < 0.15 magnitudes. 
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Table 1. Parameters of the Radial Tully-Fishcr relation at different radii for the Mathcwson sample 
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Table 2. Parameters of the Radial TF relation at different radii for the Courteau sample 
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6 APPENDIX 

In this appendix we present tables and figures related to our result and the galaxies of 
Samples 2 and 3. 
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Table 3. Parameters of the Radial Tully-Fisher relation at different radii for the Vogt sample 
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Table 4. Parameters of the standard Tully-Fisher relation for 3 samples. 
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Names of galaxies from the Courteau sample: UGC 10096, UGC 10196, UGC 10210, UGC 10224, UGC 1053, UGC 10545, 
UGC 10560, UGC 10655, UGC 10706, UGC 10721, UGC 10815, UGC 11085, UGC 11373, UGC 1152, UGC 11810, UGC 12122, 
UGC 12172, UGC 12200, UGC 12294, UGC 12296, UGC 12304, UGC 12325, UGC 12354, UGC 12598, UGC 12666, UGC 
1426, UGC 1437, UGC 1531, UGC 1536, UGC 1706, UGC 1812, UGC 195, UGC 2185, UGC 2223, UGC 2405, UGC 2628, 
UGC 3049, UGC 3103, UGC 3248, UGC 3269, UGC 3270, UGC 3291, UGC 3410, UGC 346, UGC 3652, UGC 3741, UGC 
3834, UGC 3944, UGC 4232, UGC 4299, UGC 4326, UGC 4419, UGC 4580, UGC 4779, UGC 4996, UGC 5102, UGC 540, 
UGC 562, UGC 565, UGC 5995, UGC 6544, UGC 6692, UGC 673, UGC 7082, UGC 732, UGC 7549, UGC 7749, UGC 7810, 
UGC 7823, UGC 783, UGC 784, UGC 8054, UGC 8118, UGC 8707, UGC 8749, UGC 8809, UGC 890, UGC 9019, UGC 9366, 
UGC 9479, UGC 9598, UGC 9745, UGC 9753, UGC 9866, UGC 9973 

Names of galaxies from the Vogt sample: UGC 927, UGC 944, UGC 1033, UGC 1094, UGC 1437, UGC 1456,UGC 1459, UGC 
2405, UGC 2414, UGC 2426, UGC 2518, UGC 2618, UGC 2640, UGC 2655, UGC 2659, UGC 2700, UGC 3236, UGC 3270, 
UGC 3279, UGC 3289, UGC 3291, UGC 3783, UGC 4275, UGC 4287, UGC 4324, UGC 4607, UGC 4655, UGC 4895, UGC 
4941, UGC 5166, UGC 5656, UGC 6246, UGC 6437, UGC 6551, UGC 6556, UGC 6559, UGC 6718, UGC 6911, UGC 7845, 
UGC 8004, UGC 8013, UGC 8017, UGC 8108, UGC 8118, UGC 8140, UGC 8220, UGC 8244, UGC 8460, UGC 8705, UGC 
10190, UGC 10195, UGC 10459, UGC 10469, UGC 10485, UGC 10550, UGC 10981, UGC 11455, UGC 11579, UGC 12678, 
UGC 12755, UGC 12792, UGC 150059, UGC 180598, UGC 210529, UGC 210559, UGC 210629, UGC 210634, UGC 210643, 
UGC 210789, UGC 211029, UGC 220864, UGC 221174, UGC 221206, UGC 251400, UGC 260640, UGC 260659, UGC 330781, 
UGC 330923, UGC 330925, UGC 330996, UGC 331021 
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